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МОДЕЛЬНЫЙ РАСЧЕТ ЭЛЕКТРОННО-ЭНЕРГЕТИЧЕСКОЙ 
СТРУКТУРЫ АЛМАЗА, ВМ И ВеСМ, 


Исследована электронно-энергетическая структура ряда подобных соединений: 
алмаза (С), ВМ и ВеСМ>. Установлено, что для ВеСМ› вершина валентной полосы 
формируется р-состояниями №, а дно полосы проводимости - $-состояниями Ве и С. 
Оценена ширина запрещенной полосы Е, в ВеСМ>. Проведено сравнение парциаль- 
ных и полной плотностей состояний алмаза и ВМ с рентгеновскими К-спектрами 
эмиссии и поглощения и рентгеноэлектронными спектрами, а также сравнение со- 
ответствующих парциальных плотностей состояний последних с соответствующими 
плотностями в ВесМ.. 

Ключевые слова: электронно-энергетическая структура, метод локального коге- 
рентного потенциала, плотности электронных состояний 


Введение. Поиск новых веществ со свойствами, близкими свойствам из- 
вестных материалов, а также с превосходящими параметрами был всегда 
важен с точки зрения прикладной науки. Отправной точкой поиска часто 
служат хорошо изученные и широко применяемые соединения, а новые 
получаются в результате внесения различных искажений (допирование, 
вакансии, замещение и т.д.) в их кристаллическую решетку. Теоретическое 
исследование новых (быть может, еще не синтезированных) соединений 
должно базироваться на сравнении энергетических структур этих и уже 
известных подобных соединений. 

Алмаз и ВМ уже давно представляют интерес с точки зрения при- 
кладного применения как сверхтвердые материалы, а также и с точки зре- 
ния академического интереса к их электронно-энергетическому строению 
(39С). Отметим некоторые из работ, посвященые исследованию ЭЭС. Элек- 
тронно-энергетическое строение в алмазе исследовалось методами рентге- 
ноэлектронной спектроскопии [1-3], а также рентгеновской эмиссионной 
[4] и абсорбционной [5] спектроскопии. Так, в одной из последних работ 
[3] приводится экспериментальное значение 23.0 + 0.2 эВ ширины валент- 
ной полосы в алмазе, которое сравнивается с теоретическим значением 
ширины в 21.5 + 0.2 эВ, полученным из функционала локальной плотности. 
Это различие в 2 эВ объясняется наличием возбужденной дырки и связан- 
ными с ней многоэлектронными эффектами. 

Соединению ВМ уделено в литературе и продолжает уделяться 
также достаточно большое внимание. Это относится к рентгеноспектраль- 
ным исследованиям его гексагональной модификации [6-8], а также анало- 
гичным исследованиям его кубической модификации [7, 9, 10]. Отметим 
также лишь некоторые теоретические работы по исследованию электрон- 
но-энергетической структуры ВМ [11-13]. Однако результаты работ [11, 12] 
не доведены до кривых плотностей состояний (приводятся только зависи- 
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мости Е(К) и карты зарядовых плотностей) и их трудно сравнивать с соот- 


ветствующим экспериментом по электронно-энергетическому строению. 

Несомненный интерес в последнее время появился к алмазу и ВМ и 
с точки зрения допирования их атомами с целью изменения физических 
свойств. Например, в работе [14] достаточно подробно сделано сравнение 
теоретических предсказаний уровней энергии в результате допирования с 
экспериментальными данными в ряде соединений, в том числе и для 
алмаза. 

В настоящей работе исследуется гипотетическое соединение ВеСМЬ, 
упоминание о котором было обнаружено в работе Джефа и Зангера [15], и 
его электронно-энергетическая структура (Э9С) сравнивается с ЭЭС подоб- 
ных соединений - алмаза и ВМ. По мнению авторов работы [15], соедине- 
ние ВеСМ› должно иметь свойства, приближающиеся к свойствам самого 
алмаза, а ширина запрещенной полосы Е’ должна превзойти значение Е. в 
алмазе (Е, = 5.2 эВ [16]) и быть равным 8.2 эВ. Чтобы рассчитать элек- 
тронно-энергетическое строение какого-либо соединения, необходимо чет- 
ко представлять его кристаллическую структуру. Какой она будет для гипо- 
тетического соединения можно только предполагать. Наиболее вероятной 
для ВеСМ› является структурный тип халькопирита, что может быть обос- 
новано следующими рассуждениями. Элементарную ячейку структуры ал- 
маза можно рассматривать как две гранецентрированные подрешетки, за- 
полненные атомами углерода, сдвинутые на четверть главной диагонали. 
При переходе к ВМ (изоэлектронному аналогу С) одна из подрешеток за- 
полняется атомами бора, другая - атомами азота. Получается структура 
сфалерита (цинковой обманки). Если теперь подрешетку бора представить 
как состоящую из двух регулярных подрешеток, одна из которых заполнена 
атомами Ве, а другая - атомами С, то получим структуру халькопирита. 


Пространственная группа халькопирита 1424(2’), в его примитивной 


единичной ячейке имеется 8 атомов. Кристаллическую решетку халькопи- 
рита можно считать сверхрешеткой структуры сфалерита с пространствен- 


ной группой Е43(Т?) и двумя атомами на примитивную единичную 


ячейку. В соединении ВеСМ› каждый атом азота координируется двумя 
атомами бериллия и двумя атомами углерода, а в свою очередь каждый 
атом Ве и каждый атом С координируется четырьмя атомами №. Тетраго- 
нальная элементарная ячейка халькопирита определяется тремя кристал- 
лографическими параметрами: постоянными кристаллической решетки а и 


с (чаще вместо с берется значение п = © )и параметром и, который ха- 
2а 


рактеризует смещение атома М в случае ВеСМ, из идеального тетраэдриче- 
ского положения в структуре сфалерита. Так как экспериментальные зна- 
чения а, п, и для ВеСМ, отсутствуют, то в настоящей работе для их расчета 
использовалась теория Джефа и Зангера [15]. Теория построена на двух 
предположениях: 

1) длина связи между атомами определяется как сумма их атомных 
радиусов (в данной работе использовались тетраэдрические радиусы По- 
линга [17]); 

2) предполагается, что углы связи атома углерода с четырьмя бли- 
жайшими атомами азота являются такими же, как и в идеальном тетраэдре. 
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Для ВеСМ› получены следующие значения кристаллографических 
параметров: а = 0.3847 нм; с = 0.67938 нм; п = а = 0.883; и = 0.313. 


Именно с этими параметрами производился расчет электронно- 
энергетической структуры этого соединения. 
Метод расчета. Расчеты электронно-энергетической структуры соедине- 
ния ВеСМ› проводились с использованием кластерной версии приближения 
локального когерентного потенциала [18], основанного на формализме 
функций Грина в методе многократного рассеяния. Для сравнения с ВеСМ, 
были рассчитаны также кубические модификации алмаза (пространствен- 


ная группа ЕЧЗт(О, ) ; а = 0.356688 нм [16]) и соединения ВМ (простран- 
ственная группа Е43 (Ту); а = 0.36155 нм [16]). 


Метод расчета плотностей электронных состояний (005$), предло- 
женный Джиорфи и Стотгом [18], является модификацией метода много- 
кратного рассеяния и назван методом локального когерентного потенциала 
или вторым приближением в схеме средней {-матрицы. Отличительной 
чертой метода является возможность принимать во внимание искажения 
кристаллической решетки. 

Кристаллический ти п-4п потенциал был рассчитан согласно под- 
ходу, представленному в работе Вольфа, Дякина и Широковского для кри- 
сталлической решетки с базисом [19]. Обменный член кристаллического 
потенциала рассчитывался в Х,-приближении. При этом выбирались сле- 
дующие значения о для атомных сфер различных атомов соединения и в 
области между сферами: для всех атомов и "пустых" сфер « = 1.0; для об- 
ласти между сферами «и = 0.67. 

Так как решетки сфалерита и халькопирита не являются плотно- 
упакованными, то для уменьшения объема элементарной ячейки вне атом- 
ных сфер, где потенциал имеет постоянное значение, используется стан- 
дартный прием - добавление "пустых" сфер с координатами (0.5, 0.5, 0.5) 
и (0.75, 0.75, 0.75) в единицах постоянной решетки [20]. Эти "пустые" сфе- 
ры имеют неэквивалентные положения (разные окружения) и должны рас- 
сматриваться как отдельные типы, что приводит в результате к значитель- 
ному возрастанию времени счета при большом кластере. В наших расчетах 
использовались кластеры, содержащие около 150 атомов (вместе с "пусты- 
ми" сферами). 

Для определения локальной парциальной РО$ каждого типа атома 
строился отдельный кластер. Например, в соединении ВМ для атомов В и М, 
а также пустых сфер 01 и 02, расположенных в тетраэдрических и октаэд- 
рических пустотах кристаллической решетки сфалерита, был сконструиро- 
ван соответствующий кластер, включающий примерно 150 рассеивателей. 
Центральный узел был занят (в противоположность методу средней +- 
матрицы) атомами В, М или пустыми сферами О1 и О2 соответственно. 

Плотность электронных состояний атома А-типа определяется вы- 
ражением: 


УЕ“? пиТетА® (Е) 


в (Е) =^— ] [ВА (Е) Раг А (Е) 
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где # - орбитальное квантовое число; В^(Е,г) — радиальная часть волно- 


у . А,00 м 
вой функции электрона атома А; Т,„’’ — матричный элемент опера- 
тора рассеяния; :^ — одноузельная #— матрица рассеяния для сфери- 


чески симметричных ти п-бп потенциалов; г» — радиус атомной 
сферы, равновеликой ячейке Вигнера-Зейтца. 
Полная плотность электронных состояний, приходящаяся на эле- 
ментарную ячейку, может быть найдена тогда следующим образом (напри- 
мер, для ВМ): 


(Е) = Упр (Е) + Хп^ (В) + ХВ) + вр (В). 


Сходимость метода была описана в деталях в работе [21]. Пред- 
ставленный выше расчетный метод уже применялся нами успешно к трой- 
ным халькогенидам со структурой халькопирита [22]. 

Результаты и их обсуждение. Поскольку для гипотетического 
соединения ВеСМ› нет экспериментальных данных по электронно- 
энергетической структуре, то результаты проведенных расчетов для ВеС\> 
сравниваются с выполненными в настоящей работе расчетами кривых 
плотностей электронных состояний для алмаза и кубического (структура 
сфалерита) ВМ. Сравнение рассчитанных плотностей электронных состоя- 
ний в алмазе (рис.1) и ВМ (рис.2) с экспериментальными кривыми носит в 
настоящей работе скорее тестовый характер, хотя для ВМ были взяты дос- 
таточно современные экспериментальные данные. Все рассчитанные кри- 
вые плотностей электронных состояний на рис.1,2 совмещались с соответ- 
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Рис.1. Рассчитанные плотности электронных р-состояний углерода 
и экспериментальная рентгеновская эмиссионная СК-полоса углерода [4] 
и СК-край поглощения углерода [5] в алмазе 


394 


Вестник ДГТУ, 2009. Т9. №3(42) 








ствующими экспериментальными кривыми по максимуму А, а на рис.3,4 
совмещение теоретических кривых было сделано по середине запрещенной 
полосы (провал на плотности состояний был идентифицирован с запре- 
щенной полосой Е, исследуемого соединения). Из сделанного сравнения на 
рис.1 рассчитанной плотности р-состояний углерода с рентгеновскими К- 
спектрами эмиссии [4] и поглощения [5] углерода, отражающими распре- 
деление по энергии вышеуказанных состояний соответственно в валентной 
полосе и полосе проводимости, следует, что теоретическая кривая весьма 
хорошо соответствует экспериментальным кривым. 

Для бинарного соединения ВМ плотности электронных р-состояний 
обоих компонентов соединения, т.е. бора и азота, сравниваются с соответ- 
ствующими экспериментальными спектрами из [10] на рис.2. Аналогично 
тому, что получено на рис.1 для алмаза, можно констатировать весьма хо- 
рошее соответствие расчета эксперименту. Элементы структуры А, А', В на 
теоретической кривой М р-плотности состояний совпадают по энергии с 
аналогичными элементами А, А', В экспериментальной М К-полосы. 
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Рис.2. Рассчитанные плотности электронных р-состояний бора и азота 
и экспериментальные рентгеновские К-спектры эмиссии (ВК- и МК- полосы) и поглощения 
(ВК- и МК- края поглощения) бора и азота [10] в кубическом ВМ (структура сфалерита) 


То же самое можно сказать о заполненных р-состояниях бора, для 
которых теоретический расчет воспроизводит почти всю эксперименталь- 
ную В К-полосу, за исключением максимума С. Это связано с тем, что дан- 
ный расчет в большинстве случаев не может отображать все состояния ва- 
лентной полосы и часть из них остается под ти п-йп нулем, что является, 
как известно [23], особенностью использования п"-потенциала. Весьма 
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богатая структура экспериментальных К-спектров поглощения В и М не вся 
отражена на теоретических кривых, хотя основные детали можно считать 
полученными. Получение острого пика а на ВК-крае, который авторы [10] 
связывают с рентгеновским экситоном, и вовсе невозможно в настоящем 
одноэлектронном расчете. 

Локальные парциальные плотности электронных состояний в гипо- 
тетическом халькопирите ВеСМ› приведены на рис. 3 и 4. Здесь они срав- 
ниваются с соответствующими плотностями состояний в алмазе и ВМ. Так, 
на рис.3 р-состояния углерода в ВеСМ› приведены ниже р-состояний угле- 
рода в алмазе. Происходит перераспределение интенсивности максимумов 
в валентной зоне: максимум А уменьшил свою интенсивность, а максимум В 
наоборот увеличил. Однако не наблюдается катастрофического изменения 
формы сравниваемых кривых для р-состояний углерода в алмазе и ВеСМ.. 
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Рис.3. Сравнение рассчитанных плотностей электронных р-состояний в алмазе, 
кубическом ВМ (структура сфалерита) и гипотетическом ВеСМ» 
(структура халькопирита) 


Еще меньшему изменению подверглись валентные р-состояния №, 
также приведенные на рис.3, а для свободных М р-состояний характерно 
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лишь некоторое изменение интенсивностей всех пиков. Сравнивать р- 
состояния Ве в ВеСМ› можно только с р-состояниями С и М в этом же со- 
единении, что и сделано на рис.3. Занятые р-состояния Ве лежат у верши- 
ны валентной полосы, а свободные р-состояния Ве дают интенсивный пик с 
примерно на 10 эВ выше дна зоны проводимости. В целом же можно кон- 
статировать из проведенного сравнения на рис.3, что верхнюю часть ва- 
лентной полосы в ВеСМ› формируют р-состояния М, а р-состояния Ве и С 
примешиваются к этим состояниям. Для дна зоны проводимости ВеСМ, ха- 
рактерно наличие р-состояний всех компонентов: Ве, С, М, хотя самое дно 
занимают р-состояния С, далее идут р-состояния Ве, а р-состояния М при- 
мешиваются к ним. 

Для приведенных на рис.4 кривых, отображающих $-состояния всех 
компонентов в ВеСМ›, характерно то, что занятые $-состояния Ве лежат у 
вершины валентной полосы, а аналогичные им $-состояния С (очень похо- 
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Рис.4. Сравнение рассчитанных плотностей электронных $-состояний в алмазе, 
кубическом ВМ (структура сфалерита) и гипотетическом ВеСМ» (структура халькопирита) 
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жи на $-состояния С в алмазе) и $-состояния М (не совпадающие с 5- 
состояниями М в ВМ) смещены вглубь валентной полосы примерно на 5 ЭВ. 
Свободные $-состояния С (очень похожие по форме на $-состояния С в ал- 
мазе) и $-состояния Ве в ВеСМ› расположены у дна зоны проводимости. На 
это указывает энергетическое положение пиков а и Ь на кривых С 5- 
состояний и Ве $-состояний в ВеСМ.. На самом дне зоны проводимости на- 
ходятся $-состояния М. В целом по результатам расчета можно говорить о 
высокой степени гибридизации как свободных р- и $-состояний, так и заня- 
тых р- и $-состояний компонентов в ВесСМ.. 

Весьма важным, на наш взгляд, является сравнение кривых, приве- 
денных на рис.5, где показаны полные плотности состояний в алмазе, ВМ и 
ВеСМ>. Как видно, валентная полоса в ВеСМ› шире, чем в ВМ. Кривая сво- 
бодных состояний в ВеСМ› имеет ряд элементов структуры, которые сбли- 
жают ее с формой кривой свободных состояний в алмазе (элементы Ь и с) 
ив ВМ (элемент с"). Во всех трех случаях (см.рис.5) имеется область энер- 
гий, где плотность спадает до нуля, и эта область, несомненно, представля- 
ет оптическую запрещенную полосу Е. в этих полупроводниках. Для алмаза 
и ВМ это подтверждается сравнением с экспериментальными данными на 
рис.1 и 2, где рентгеновские эмиссионные и абсорбционные спектры доста- 
точно хорошо совпали по основным элементам структуры с рассчитанными 
здесь теоретическими плотностями состояний и тем самым позволили 
четко определить область запрещенных энергий или так называемую 
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Рис.5. Рассчитанные полные плотности электронных состояний в алмазе, 
кубическом ВМ (структура сфалерита) и гипотетическом ВеСМ› 
(структура халькопирита) 
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оптическую запрещенную полосу Е’ на теоретической кривой. Оценка ши- 
рины запрещенной полосы ЛЕ проведена для всех исследованных полупро- 
водников и сравнивается с литературными Е; в таблице. 


Кристаллографические параметры, ширины оптических запрещенных полос 
Е, и оценочные значения ДЕ для исследованных полупроводников 





























а, нм С, НМ Е, ЭВ ДЕ, эВ 

Алмаз 0.356688 [16] : 5.2 [16] 6.0 +0.5 
- 6.00.5 [16] 

вм 0.36155 [16] а 6505 

Вес\ 0.3847 0.67938 к . 





Оценочные значения получились несколько завышенными по срав- 
нению с литературными, хотя сами литературные данные имеют значи- 
тельный разброс. Отметим, что полученное нами значение Е. = 6.0 + 0.5 эВ 
в ВеСМ› отличается от предсказанного в работе [15] значения в 8.2 эВ на 
2.2 эВ, что уже нельзя просто связать с неточностью метода определения. 
Выводы. Проведенные расчеты показывают качественное подобие элек- 
тронно-энергетической структуры (3ЭС) гипотетического тройного 
соединения ВеСМ) его изоэлектронным и изоструктурным аналогам — алма- 
зу и ВМ. 

Основные особенности кривых плотности электронных р-состояний 
М в ВМ и в ВеСМ> достаточно хорошо соответствуют друг другу. Форма ва- 
лентных р-состояний М в ВеСМ, практически не отличается от формы ана- 
логичных состояний в ВМ. В то же время форма кривой в области свобод- 
ных р-состояний № в ВеСМ› претерпевает существенные изменения по 
сравнению с ВМ. 

В отличие от р-состояний кривые плотности $-состояний М в ВеСМ» 
претерпевают качественные изменения по отношению к ВМ. Происходит 
перераспределение плотности как занятых, так и свободных $-состояний № 
в ВеСМ› по сравнению с ВМ. Подобным же образом наблюдаются качест- 
венные отличия в форме кривых р-состояний С в алмазе и ВеСМЬ, тогда как 
форма кривых их $-состояний очень похожа. 

Расчеты показали подобие формы кривых р-состояний Ве в ВеСМЬ 
р-состояниям С в алмазе и подобие формы кривых 5-состояний Ве 5- 
состояниям М в кубическом ВМ. 
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А.А. ГАУВЕМТУЕМ, В.У. САВКЕГТАМ, Т.Уа. МТКТРОВОМУ 


МОБЕЕ САЕСУТАТТОМ ОЕТНЕ ЕТЕСТКОМ ЕМЕВСУ 5ТОСТОВЕ 
ОЕ ОТАМОМО, ВМ А\О ВеСМ, 


Тре весгопс епегду гисиге оР {Пе го\м/ оЁ этйаг сотрочуп@$: Фатопа (С), 
ВМ, апа ВеСМ>› Ваз Бееп туезНдакеа. Г мгаз Гоипа Гог ВеСМ» {ПаЁ Пе юр оЁ 
{пе уаепсе Бапа 15 огтеа Бу р-$ае$ оЁ М, мпегеа$ Не БоКот ог {Пе сопдис- 
Чоп Бапа - Бу 5-5ае$ оЁ Ве апа С. Тпе миа оЁ {Пе епегду ГогГладеп Бапа Ех 
ш ВесСМ› — маз езбтаеа. ТНе сотрамзоп оЁ {Ле сасШаей рафа! ап {ока! 
еесгоп депзШез оЁ $а(ез ог фатопа апа ВМ МИ х-гау К-зресга оЁ епт5$!оп 
ап абзогриоп апа мА х-гау @есгоп зресга \маз сагйе ои{. ТНе согге- 
зропатд раг4а! еесгоп ${аез ог ЧатопЯ апа ВМ меге а[5о сотраге у/И 
соггезропата ЧепзШез оЁ $(афез оЁ Вес. 
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